An Algorithm for the Solution of Impurity Diffusion under Finite Reaction Rates
نویسندگان
چکیده
An algorithm allowing a fast solution of stationary and time dependent 1-d diffusion problems for the case of finite ionization and recombination rates is presented. Results of numerical computations and an analytic stability analysis show its unconditional stability. The computational effort involved rises only linearly with the number of ionization stages included. An application to the diffusion of oxygen and iron impurities in JET is described.
منابع مشابه
A numerical investigation of a reaction-diffusion equation arises from an ecological phenomenon
This paper deals with the numerical solution of a class of reaction diffusion equations arises from ecological phenomena. When two species are introduced into unoccupied habitat, they can spread across the environment as two travelling waves with the wave of the faster reproducer moving ahead of the slower.The mathematical modelling of invasions of species in more complex settings that include ...
متن کاملPositivity-preserving nonstandard finite difference Schemes for simulation of advection-diffusion reaction equations
Systems in which reaction terms are coupled to diffusion and advection transports arise in a wide range of chemical engineering applications, physics, biology and environmental. In these cases, the components of the unknown can denote concentrations or population sizes which represent quantities and they need to remain positive. Classical finite difference schemes may produce numerical drawback...
متن کاملNumerical Study on the Reaction Cum Diffusion Process in a Spherical Biocatalyst
In chemical engineering, several processes are represented by singular boundary value problems. In general, classical numerical methods fail to produce good approximations for the singular boundary value problems. In this paper, Chebyshev finite difference (ChFD) method and DTM-Pad´e method, which is a combination of differential transform method (DTM) and Pad´e approximant, are applied for sol...
متن کاملEquidistribution grids for two-parameter convection–diffusion boundary-value problems
In this article, we propose an adaptive grid based on mesh equidistribution principle for two-parameter convection-diffusion boundary value problems with continuous and discontinuous data. A numerical algorithm based on an upwind finite difference operator and an appropriate adaptive grid is constructed. Truncation errors are derived for both continuous and discontinuous problems. Parameter uni...
متن کاملAn hp finite element method for singularly perturbed systems of reaction- diffusion equations
We consider the approximation of a coupled system of two singularly perturbed reaction-diffusion equations by the finite element method. The solution to such problems contains boundary layers which overlap and interact, and the numerical approximation must take this into account in order for the resulting scheme to converge uniformly with respect to the singular perturbation parameters. We pres...
متن کامل